Quarkonium Spectral Functions 



Agnes Mocsy 

Department of Mathematics and Science, Pratt Institute, Brooklyn, NY 11205, USA 



Abstract 

In this talk I summarize the progress achieved in recent years on the understanding of quarko- 
nium properties at finite temperature. Theoretical studies from potential models, lattice QCD, 
and effective field theories are discussed. I also highlight a bridge from spectral functions to 
experiment. 



1. Introduction 

It has been 23 years since what is today the most cited paper in heavy ion physics [ 1 ] was pub- 
lised. In this work Matsui and Satz conclude that "J/if/ suppression in nuclear collisions should 
provide an unambiguous signature of quark-gluon plasma formation". This is how the story 
of quarkonium at finite temperature began. One of the important and necessary consequences 
of this story is the quest for determining properties of quarkonium states at finite temperature. 
Quarkonium properties can be conveniently studied through spectral functions. In this talk I dis- 
cuss properties of quarkonium inside an equilibrated plasma, as obtained from the three main 
theoretical approaches for determining the spectral functions. These are, in historic order, poten- 
tial models, lattice QCD, and effective field theories. I then discuss what possible implications 
for the experiments the results might have. 

Quarkonium states, made of a heavy quark and its antiquark (niQ » Aqcd), are tightly 
bound; their size is much smaller than the typical hadronic scale, rj^ =* 0.4 fm and ry =* 0.2 fm 
<K 1 fm. Due to the heavy quark mass the quark velocity is very small, v <k 1, which allows 
for a non-relativistic treatment of the bound states. Properties of quarkonium have been custom- 
arily obtained by solving the Schrodinger equation, where the interaction between the quark Q 
and antiquark Q at distance r away from each other is implemented as a potential, V(r). This 
prescription works very well at zero temperature and it has been somewhat naively used also at 
finite temperature. There, the idea is that in a high temperature deconfined medium there is a re- 
arrangement of color around the heavy quarks, due to which the effective charge of the Q and Q 
is reduced. So light quarks and gluons screen the potential between the Q and the Q, analogous 
to the known Debye-screening in QED plasmas. The original Matsui-Satz idea is that when the 
range of screeing becomes smaller than the size of the state, rp < r Q -Q, the quark and antiquark 
cannot see each other and the two are no longer bound. Customarily, these effects were described 
with a temperature-dependent potential V(r, T). 

This simple and attractive idea leaves a few questions: How good is it to describe medium 
effects with a temperature dependent potential? Is there Debye-screening in the plasma? Is 
screening the only relevant effect? And what do we mean by a bound state at finite temperature? 
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2. Spectral Functions 



Quarkonium properties are conveniently described using spectral functions. These contain 
all the relevant information about a given channel, providing a unified treatment of bound-states, 
continuum, and threshold. Bound states show up as peaks in the spectral function. At zero tem- 
perature states are well defined, the widths of the peaks can be narrow, the lifetimes of the states 
can be large. The energy gap between the continuum threshold and a peak position defines the 
binding energy, Ebi„ for that state. Peaks are expected to broaden with increasing temperature. 
The disappearance of a peak from the spectral function means the bound state is dissociated (of- 
ten referred to as the state is melted). Here I want to emphasize, that bound states may disappear 
well before the binding energy vanishes. Thus a melting (or dissociation) temperature should be 
defined by a less severe condition than zero binding energy, contrary to what has generally been 
considered in the literature. 

There are three main theoretical approaches to determine spectral functions. 1) Finite tem- 
perature potential models date back to the Matsui and Satz paper. Besides they are easy to 
handle, a great advantage of the newer versions of potential models is that they easily acco- 
modate non-perturbative information from lattice QCD calculations of the free energy of static 
quarks, in form of different "lattice-based potentials". The drawback is, however, their ad-hoc 
nature. 2) The first lattice QCD calculations for quarkonium spectral functions at finite tempera- 
ture appeared in 2003. Although the computation of quarkonium correlators is exact, it is really 
difficult to extract spectral functions on quantitative and sometimes even qualitative level from 
these. 3) Effective field theories have only appeared in the last few years. These are obtained 
directly from QCD and they can clarify the domain of applicability for the simpler potential 
model approach. However, currently derivations of the potentials from effective theories must 
assume some scale separation (between the bound state- and temperature-scales), and utilize 
weak-coupling techniques, which makes their application near the critical temperature problem- 
atic. With the advantages and draw-backs of each of these approaches, a combination of all their 
results provides a newly emerging complete picture of quarkonium at finite temperature. The 
chart on figure[T]illustrates the interconnectedness of these approaches. 

3. Color Screening from Lattice QCD 

A rigorous way to study the modification of inter-quark forces with temperature T and quark 
separation r is provided from lattice QCD by computing the heavy quark singlet free-energy, 
F\(r, T) [2]. This is shown in the left panel Figure[2] At short distances, r < r me d(T), the singlet 
free energy is independent of temperature and coincides with the zero temperature potential. The 
range r me d(T) where F\{r, T) = F\(r) decreases with increasing temperature. At large distances 
r > r scr (T) and temperatures above the critical T c the singlet free-energy is screened, F\(r, T) = 
Fi(T). With increasing T the strong modification of the static Q - Q interaction sets in at shoter 
and shorter distances r scr (T). As illustrated on the right panel of Figure [2] the range of Q - Q 
interaction is comparable to the size of QQ charmonium states. This generic observation together 
with the original Matsui-Satz argument suggests that the charmonium states, such as the Jj^i, 
and excited bottomonium states, do not exist as bound states above T c , while the ground state 
bottomonium, T, is small enough to survive deconfinement up to higher temperatures. One 
problem with this line of reasoning is that the singlet free energy, Fi(r,T), in general, does 
not coincide with the potential V(r, T). It has thus been argued to consider the internal energy, 
obtained from the free energy, as the heavy quark potential. However, the internal energy is also 
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Figure 1: Flow chart of in-mdeium quarkonium calculations. 



not the potential. The right panel of Figure [2] further illustrates that at a given temperature the 
internal energy is more binding than the corresponding free energy, leading possibly to survival 
of quarkonium states when identified as the heavy quark potential. In what follows, I will discuss 
spectral functions obtained from a variety of lattice-based potentials, and argue that survival of 
the J/ifj well above T c is unlikely. 



4. Spectral Functions from Potential Models 

In potential models the basic assumption is that all medium effects are accounted for as a 
temperature-dependent potential. The potential that describes heavy quark-antiquark interaction 
in a hot medium is not yet known. Furthermore, even the applicability domain of the traditional 
potential model approach is questionable. The newest advancements towards deriving the po- 
tential from QCD are from effective field theories, discussed in a following section. Lacking an 
exact form, potentials are either phenomenological ^ or lattice-based 0. El H S H & O • 

The main advantage of lattice-based potentials is that non-perturbative effects can be conve- 
niently accommodated using lattice data on the heavy quark-antiquark singlet free energy dis- 
cussed above and shown in the left panel of Fig. [2] The true potential at a given temperature and 
large distances (if it can be addressed at all in this approach) lies somewhere between the free- 
and the internal-energy (see right panel of Fig. [2j. As said in the previous section, at short dis- 
tances it corresponds to the zero temperature potential, but we do not have information about its 
behavior at intermediate distances. Utilizing these constraints one can interpolate at intermediate 
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Figure 2: Heavy quark singlet free energy versus quark separation calculated in 2+1 flavor QCD for different tempera- 
tures 1 2] (left) and constraints from lattice on the potential at a given temperature (right). 



distances and construct lattice-based potentials. Clearly, the shortcoming of this approach is its 
somewhat ad hoc nature. We can however, look at the most binding potential still consistent with 
lattice and learn about the quarkonium spectral functions in this upper limit. 

Figure[3]shows quarkonium spectral functions obtained from a non-relativistic Greens func- 
tion calculation using the most confining potential in full QCD [10] (left panel) and from a T- 
matrix calculation using the internal energy as potential in quenched QCD [|8|] (right panel). Both 
results are for the S-wave charmonium channel. The spectral function in the left panel shows no 



CF((Ju)/(JU 2 12 
1 

0.3 

0:6 

0.4 
0.2 




Mis) 



n'c(2Sj 




' T=0 ' 
I.ITc 



I.STc 
free 



continuum 



3.5 4 4.5 

co[GeV] 




Figure 3: S-wave charmonium spectral functions calculated in potential model from fTofl (left) and from fl] (right). 



peak structures at high temperature, indicating that all charmonium states are dissolved in the 
deconfined phase JiTOH. The result shown in the right panel finds the ground state peak at 
somewhat higher temperatures yj]. 

Even though there are some differences in the details, there are some essential features com- 
mon to all models: For one, there is a large threshold (rescattering) enhancement beyond what 
corresponds to free quark propagation. This enhancement is present even at high temperatures 
and it indicates correlations persisting between the quark and antiquark even in the absence of 
bound states. Threshold enhancement has been identified in all of the channels (charmonium and 
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bottomonium S- and P-states) J9L fToll . Second, there is a strong decrease with increasing tem- 
perature of the binding energies (the distance between peak position and continuum threshold 
in the spectral function) determined from potential models [11]. The decreased binding energy 
implies a large increase in the thermal width as well. As mentioned earlier in this talk, bound 
states disappear before the binding energy goes to zero Hill . One can think that this happens 
when the width becomes larger than the binding energy. In other words the time it takes for a 
quark and antiquark to bind is larger than the time it takes for the bound state to decay. Figure H] 
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Figure 4: S-wave bottomonium (left) and charmonium (right) spectral functions calculated in potential model |9]. Insets: 
correlators compared to lattice data. Charmonium is compared to lattice data from 1 20] 



shows a ground state bottomonium peak up to higher temperatures into the deconfined medium, 
but with a largely reduced binding energy. I note that binding energies are somewhat increased in 
a viscous medium [12]. For further details on anisotropy effects on quarkonium I refer the reader 
to lll3ll . I further note that potential model calculations of spectral functions do not include a re- 
alistic temperature-dependent width for the states. Widths naturally arise from effective theories 
lHHEHL as discussed in a later section. We calculated the thermal widths from the binding ener- 
gies 11 1 Ofl following [1611. Thermal broadening of quarkonium has been addressed also in a NLO 
perturbative QCD II 1 711 and a QCD sum rule II 1 811 analysis. All of these calculations show that the 
J I ^ for example, is significantly broadened at temperatures right above that of deconfinement. 

From the analysis of the spectral functions one can provide an upper bound for the disso- 
ciation temperatures, i.e. the temperatures above which no bound states peaks can be seen in 
the spectral function. Above these temperatures bound state formation is suppressed. Conserva- 
tive upper limit dissociation temperatures for the different quarkonium states have been obtained 
from a full QCD calculation If 1 Ofl . Accordingly, charmonium states are dissolved in the decon- 
fined phase, while the bottomonium ground state may persist up to temperature of about 2T C . 



5. Spectral Functions from Lattice QCD 



In lattice QCD, current-current correlation functions of mesonic currents in Euclidean-time 
are calculated [19, 2(1 21]. These correlators are related to the spectral function through the 
integral: 



G{r. 



,71 = J 



da>o-(a>, r)cosh(w(r- 1/2T)) / sinh(<y/2r) . 



(1) 
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It is customary to present the temperature-dependence of the correlators as the ratio G/G recon , 
where G recon (T, T = 0) is the correlator reconstructed from spectral function at zero (or some low) 
temperature. Since correlators are computed with high accuracy, and they can be determined in 
potential models by integrating the spectral function using equation it has been suggested 
to compare correlators from potential models to the ones from lattice QCD [5fl- Somewhat sur- 
prisingly, Euclidean correlation functions show very little temperature dependence, irrespective 
of whether a state is there (such as the T) or not (such as the J/ iff). Note also, that correlators 
from potential models are in accordance with the lattice calculations (see insets in Fig. [3) (for a 
review see 111 111 ). Originally, the small temperature dependence of the correlators was considered 
as evidence for survival of different quarkonium states [T^, 22[] . It is now clear that this conclu- 
sion was premature. So flat correlators do not tell about the survival or melting of states, since 
these kind of changes do not show up in it. We now understand that the threshold enhancement 
compensates for the absence of bound states and leads to Euclidean correlation functions which 
show only very small temperature dependence ^ . 

From the current-current correlation functions for different quarkonium channels the quarko- 
nium spectral functions are extracted $19 , 22 , 2(3] • The extraction of spectral function using the 
Maximum Entropy Method is still difficult due to discretization effects, statistical errors, default 
model dependence (for a review see[23]). So the uncertainties in the spectral function are signif- 
icant and details of this cannot be resolved. Moreover, a seemingly surviving ground state peak 
is in perfect agreement with a mere threshold enhancement obtained from potential models, as 
shown juxtaposed on the right panel of figure [4] Again, a 1 GeV wide peak obtained with huge 
uncertainties from lattice QCD, and peak, whose details cannot be resolved, is most likely not a 
true bound state peak. Improved calculations are on the way 112411 . 



6. Spectral Functions from Effective Field Theories 

In effective field theories the quarkonium potential is derived from the QCD Lagrangian. The 
basis of these theories is the existence of scales related to the bound state and scales related to 
the temperature. At T=0 for heavy quarks the existence of different energy scales related to the 
heavy quark mass m, the inverse size mv ~ 1 fr , and the binding energy mv 2 ~ E\, m makes 
it possible to construct a sequence of effective theories. The effective theory which emerges 
after integrating out the scale m and mv 2 is pNRQCD, which delivers the potential model at 
T — [25]. Recently this approach has been extended to finite temperature, where the presence 
of the scales T, the Debye mass mo ~ gT, and the magnetic scale g 2 T, makes computations 
more difficult. With the assumption of m » T and weak coupling g <k 1 these scales are well 
separated. Depending on how these thermal scales compare with the bound state scales, the 
different hierarchies allow for the derivation of different effective theories for quarkonium bound 
states at finite temperature [15]. A finding common to all of the theories is that the quarkonium 
potential has not just a real part but also an imaginary part. The temperature dependent imaginary 
part determines the thermal width, and thus the smearing out, of the quarkonium states in the 
spectral function. Several physical processes contribute to this thermal width. For instance the 
scattering of particles in the medium with gluons (Landau damping) [ 14], or the thermal break- 
up of a color singlet Q - Q into a color octet state and gluons (octet transition) 11511 . Smearing 
of states, leading to their dissolution, even before the onset of the exponential Debye-screening 
of the real part of the potential (see e.g. discussion in UH). 

Another discovery from the new effective field theory calculations is that finite temperature 
effects can be other than the originally thought exponential screening of the real part of the 
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potential; in the weak coupling approach thermal corrections to the potential are obtained only 
when the temperature is larger than the binding energy II 1 511 . 

The current shortcomings of effective field theories are that near T c the applicability of weak 
coupling techniques is problematic and theories must incorporate also non-perturbative effects. 



7. Relevance for experiments 



Knowing the quarkonium spectral functions in equilibrium QCD is necessary, but alone is 
not sufficient to predict their production in heavy ion collisions. In principle, there is a simple 
relation between quarkonium spectral functions and quarkonium production rates measurable in 
experiments. This relation will hold only in thermal equilibrium, which is not likely the achieved 
for heavy quarks in heavy ion collisions. The bridge between theoretical spectral functions and 
experimental yield measurements requires additional dynamical mode ling (se e figure [T]). There 
are several such models with different underlying assumptions 

Let us look for instance at RHIC energies, where we learned from potential models that no 
J/ip peak in the spectral function is seen, only threshold enhancement This means that 

J /if/ are not formed, only spatially correlated c and c JSjllOl]- This correlation can be modeled 
classically, using Langevin dynamics which includes a drag force and a random force on the 
heavy quark (antiquark) from the medium, as well as the forces acting between the quark and 
anti-quark (described by the potential). The distribution of the separation between c — c pairs 
with and without threshold enhancement in the spectral function is shown in the left panel of 
Figure |5] Some of these c — c pairs will stay correlated throughout the evolution of the system. 
Those pairs that are not diffused away they can bind together into a J /if/ at hadronization. 

The right panel of Figure [5] shows the result of a Langevin simulation of evolving c—c pairs 
on top of a hydrodynamically expanding quark-gluon plasma, which describes the RHIC data 
on charmonium suppression quite well B29I1 . In particular, this model can explain why, despite 
the fact that a deconfined medium is created at RHIC, we see only 40 - 50% suppression in the 
charmonium yield. 
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Figure 5: Distribution of c - c pairs versus the inter-quark distance (left) and predicted suppression (right) from [ 29]. 
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8. Summary 



Potential model calculations based on lattice QCD, as well as resummed perturbative QCD 
calculations indicate that all charmonium states and excited bottomonium states dissolve in the 
deconfined medium. Lattice data is consistent with J /if/ screened just above T c . Dissolved states 
lead to the reduction of the quarkonium production yield in heavy ion collisions compared to the 
binary-scaled proton-proton collisions. Due to possible recombination effects, however, the yield 
will not be zero. Potential models indicate the decrease in the binding energy of quarkonium 
states and a sizable threshold enhancement leading to residual quark-antiquark correlations per- 
sisting at high temperatures. Implication of this for the understanding of experimental data can 
be investigated in dynamical models, which suggest that recombination of residually correlated 
pairs is possible when the temperature cools down sufficiently. Effective field theory calculations 
indicate temperature-dependent widths and temperature effects beyond what potential models 
have accounted for. 
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